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Abstract 

We derive explicit reconstruction formulas for the attenuated geodesic X-ray transform 
over functions and, in the case of non-vanishing attenuation, vector fields, on a class of 
simple Riemannian surfaces with boundary. These formulas partly rely on new explicit 
approaches to construct continuous right-inverses for backprojection operators (and, in turn, 
holomorphic integrating factors), which were previously unavailable in a systematic form. 
The reconstruction of functions is presented in two ways, the latter one being motivated 
by numerical considerations and successfully implemented at the end. Constructing the 
right-inverses mentioned require that certain Fredholm equations, first appearing in [29], be 
invertible. Whether this last condition reduces the applicability of the overall approach to a 
strict subset of simple surfaces remains open at present. 


1 Introduction 

We consider explicit inversion formulas for the two-dimensional attenuated X-ray (or, equiv¬ 
alently in two dimensions, Radon) transform of a function or vector field. Let (M, g) a non¬ 
trapping Riemannian surface-with-boundary with unit circle bundle 

5M = {(x,n) gTM, g{v,v) = l}, 

and let us denote the geodesic flow tpt : SM — )• SM by tpt{x,v) = where ^x,v{t) 

denotes the basepoint and ^x,v{t) is the (unit) velocity vector. For x G dM, let Vx the unit inner 
normal and define the influx/outflux boundaries d±SM = {(x,u) G dSM, ±g{ux,v) > 0}. Fix 
a a smooth enough function on M. Then for a function / G Lp‘{SM), we dehne the attenuated 
(geodesic) X-ray transform of / by 

rT(x,v) / r-t \ 

Iaf{x,v):=J f{ipt{-K,v))ex.p Ij a(jx,v(s)) dsj dt, {x,v) e d+SM, 
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where r(x, denotes the first exit time of the geodesic 7x,t;(i) (the non-trapping condition 
implies that r is uniformly bounded above on SM). Such a transform can be regarded as the 
influx trace laf = v\d+SM of the solution u to the transport equation on SM 

Xu + au = -f (SM), u\3_sM = 0, (1) 

and where X = ^(pt\t=o is the geodesic vector field. Cases of interest here are (i) the case where 
/ is a function on M i.e. f{x,v) = /(x) with applications to X-ray Tomography in media with 
variable refractive index, a topic receiving much interest at the moment [22, 15], and {ii) the 
case where / represents a vector field F on M via the relation f{x,v) = g{F{x),v), and whose 
applications to Doppler tomography justify its nickname of Doppler transform. 

Such a transform over functions was studied earlier in the Euclidean setting. Inversion meth¬ 
ods with known attenuation were obtained independently by Arbuzov, Bukhgeim and Kazantzev 
using A-analytic function theory d la Bukhgeim in 1998 [1], and by Novikov via complexifica- 
tion methods in 2002 [25, 23], see [4] for a joint study of both approaches. While both methods 
present two interesting ways of looking at this problem, the latter is of lesser complexity (i.e., 
comparable to that of an inverse Radon transform) and is therefore more amenable to imple¬ 
mentation, as illustrated in [21, 7]. This latter method was extended to partial data and more 
general integrands in [2], to the attenuated Radon transform in hyperbolic geometry and to the 
horocyclic transform in [3], to even more general curves in [10], and more recently, to weighted 
Radon transforms in the Euclidean plane, where the weight has finite harmonic content in [24], 
with an implementation in [8]. Both methods were also adapted to the study of the attenuated 
transform over functions and vector fields in fan-beam geometry in [12], describing both a so¬ 
lution via A-analytic functions and an approach using fiberwise holomorphic solutions to the 
transport equation (1) in the Euclidean setting. This last approach was then generalized by Salo 
and Uhlmann to simple Riemannian surfaces (i.e. surfaces with strictly convex boundary and 
no conjugate points) in [34], where a method reconstructing functions from knowledge of their 
geodesic X-Ray transform was developed. The Doppler transform was studied microlocally in 
the Riemannian case in [11], and a range characterization in the Euclidean case was recently 
given in [32]. Additional range characterizations of the Euclidean transform in convex domains 
of were provided in [33] , and a study of the attenuated Euclidean transform over second-order 
tensors was recently provided in [31]. 

The X-ray transform can also be considered over vector-valued unknowns defined on SM (i.e. 
sections of bundles), with matrix-valued attenuations, connections and Higgs fields, for which 
recent results can be found in [27, 26, 6]. More general settings include the case of weighted 
X-ray transforms (of which the attenuated case is a particular example). In this setting, early 
Euclidean implementations have been provided in [14], and analytic microlocal analysis has 
led in [5] to injectivity and stability of the attenuated transform over functions when both the 
metric and the attenuation coefficient are real-analytic. In dimension d > 3, local injectivity 
of weighted X-ray transforms near convex boundary points was recently established in [37], 
following methods in [36] where the unattenuated case was first treated. 
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While laying the groundwork of inversions on surfaces with non-constant curvature in [34] by 
proving injectivity and providing the first reconstruction procedure, the authors there pointed 
out two open questions: (i) how to explicitely invert the unattenuated X-ray transform over 
functions and solenoidal vector fields (we call them Iq and I± here), and (ii) how to explicitely 
construct (fiber-)holomorphic integrating factors when curvature is not constant. Here by holo- 
morphic integrating factor (HIF) for the attenuation o, we mean a function which is holomorphic 
on the fibers (see Sec. 2.1), and which turns the attenuated transport equation (1) into an unat¬ 
tenuated one. As will be seen below, after performing Fourier analysis on the fibers of SM, a 
transport equation such as (1) takes the form of a tridiagonal, doubly infinite system of partial 
differential equations on M, which is best understood when it can be made one-sided. As fiber- 
holomorphic functions have one-sided harmonic content, HIFs allow to fulfill this requirement to 
a certain extent. In all past aproaches analyzing attenuated transforms, a “holomorphization” 
step of integrating factors always occurs, be it on the fibers of SM when solving a transport 
equation [12, 34], in terms of a complexified parameter when formulating a Riemann-Hilbert 
problem [25, 3, 2], or when considering sequence-valued systems [1, 31, 32, 33]. 

In an attempt to address the open questions mentioned above, an answer to (i) was presented 
by the author in [18], by implementing the reconstruction formulas proposed in [29, 13, 28]. A 
first feature of the present article is to address point (ii) by providing explicit ways of construct¬ 
ing holomorphic integrating factors, obtained by first constructing preimages of the adjoint 
operators Jq and I'f explicitely. Second, we derive other reconstruction algorithms, some of 
which are similar in spirit to those in [12], another one similar to the approach in [29, 19]. In 
the first approach, the main novelty here, partly motivated by the approach in [12], consists in 
decomposing v = ue~'^ (with u the solution of (1) and w a holomorphic solution of Xw = —a) 
into the sum of a holomorphic function and a function that is constant along geodesics. This can 
be done as soon as one can construct explicit preimages of the operators Iq and If. A second 
approach reconstructing functions, more similar to approaches in [29, 19], is then presented, and 
a numerical implementation is provided. Additionally, these reconstruction formulas are fast 
in that no full three-dimensional transport equation needs to be solved unlike [34]. The class 
of surfaces where the current approach is valid is that of simple surfaces where some Fredholm 
equations (see Equations (10) and (11) below), which first appeared in [29, Theorem 5.4], are 
invertible. It remains open at present whether this latter additional requirement applies to all 
or a strict subset of simple surfaces, though past numerical experiments done in [18] by the 
author have showed that such Fredholm equations were invertible on some family of surfaces 
which could become arbitrarily close to non-simple. 

Recent work by the author with P. Stefanov and G. Uhlmann [20] shows that stable inversion 
of the attenuated ray transform should still be possible in some surfaces where conjugate points 
occur at most in pairs. Adapting the current approach to this latter setting will be the object 
of future work. 
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Outline. The structure of the paper is as follows. We first introduce in Sec. 2 the basic 
setting (Sec. 2.1), some new notation and operators (e.g., I±) which play an important role 
in the inversion process (Sec. 2.2), the construction of explicit, continuous right-inverses for 
Iq and (Sec. 2.3) and how they allow the explicit construction of holomorphic integrating 
factors (Sec. 2.4). Section 3 presents the reconstruction formulas for functions and vector fields, 
following initial ideas in [12], first adapted in [34]. Section 4 proposes an alternative approach 
to reconstruction of functions, stating in passing —)• continuity and bounds for a certain 
family of operators, the proof of which is relegated to appendix A. Section 5 presents a numerical 
implementation of inversions as set up in Theorem 4.2. 

2 Notation, preliminaries and the unattenuated case 

2.1 Geometry of the unit circle bundle 

We briefly recall the geometry and notation associated with the unit circle bundle. With (M, g) 
as above, the geodesic flow (ft : SM —)• SM is defined on the domain 

= {(x,u,t) : (x,u) G SM, t G (-r(x,-u), r(x, u))}. (2) 

with X the geodesic vector field as in the introduction, one may construct a global frame 
{X,X±, V} of T{SM), which encodes the geometry of M via the structure equations 

[X, V] = X±, [A"_l, V] = X, [X, X±] = —kV {k : Gaussian curvature). 

V is sometimes referred to as the vertical derivative and X_\_ the horizontal derivative. For 
X G M, we also denote 5x := {v G T^M, g{v,v) = 1} the unit tangent circle at x. 

Though the proofs will be coordinate-free, a convenient set of coordinates is that of isothermal 
coordinates (they can be made global as the simplicity assumption on M implies that it is simply 
connected), for which the metric g is scalar of the form g = -|- dy^), and where 

(x, v) G SM is parameterized by (x, y, 0) where x = (x, y) and v = (cos 6dx + sin 9dy) for 

0 G In these coordinates, the frame {X,X±, G} takes the expression 

X = e“^(cos ddx + sin 9dy + (— sin 9dxX + cos 9dyX) dg), V = dg, 

X± = —e~^{— sm9dx + cos9dy — {cos 9dxX + sm9dyX) dg). 

We use the Sasaki metric on SM for which the basis {X, X±,V) is orthonormal, with volume 
form preserved by the frame (in isothermal coordinates, = e^^ dx dy d9). Introducing 
the inner product 


{u,v) 


uv dS^, 


Ism 


u, V : SM — >■ C, 
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the space L‘^{SM,C) decomposes orthogonally as a direct sum 

L^{SM,C) = Hk, where := kev{V — ikid). 

kez 

We also denote := C°°{SM) n and denote the corresponding decomposition u{x,v) = 
with Uk e Hk- In isothermal coordinate, this corresponds to Fourier series expan¬ 
sions 

“ If 

u{x,9)= Uk{x,0), where nfc(x, 0) = nfc(x) = — / u{x,9)e~^^^d9. 

k=—oo 

In particular, we denote either by uq or vrou the hberwise average no(x) = ^ Jg u(x,v) dS{v), 
so vTo : L?‘{SM) —)• L?‘{M) is the projection onto Hq. Such functions admit an even/odd 
decomposition w.r.t. to the involution v i—>■ —v (or 9 ^ 9 + n), denoted 

u = + U-, where := Uk and U- : = E Uk- (4) 

k even k odd 

An important decomposition of X and X± due to Guillemin and Kazhdan (see [9]) is given by 
defining rj± := so that one has the following decomposition 

X = r]^+ri_ and X_l = t (r?-i-- ??-), 

i 

with the important property that r]± : Q,k ^ Hfeii for any k € Z, so that both X and X± map 
odd functions on SM into even ones and vice-versa. 

In the harmonic decomposition above, a diagonal operator of particular interest is the so- 
called fiberwise Hilbert transform H : C°°{SM) —)• C°°{SM), whose action on each component 
is described by 

Huk '■= —isign{k)uk, k £ Z, with the convention sign(O) = 0, 

and we denote H_^/_ the composition of H with projection onto even/odd Fourier modes. We 
say that a function u G L‘^{SM) is (fiber-Jholomorphic if {Id + iH)u = uq, i.e. if u has only 
nonnegative Fourier components. An important identity first proved in [30] is the commutator 
between the Hilbert transform and the geodesic vector field: 

[H, X] = ttoAx + [H, Ax] = -ttoA - Atto. 

Note also that = —Id -\- vro and that Ht^q = ttqH = 0. Using these observations and the 
commutators above, we write 

[if^ A] = H[H, A] + [H, X]H = IHftr^+ HX^tto + + ttqX^H 

On the other hand, [H‘^,X] = [—Id vro. A] = vroA — Avro. Upon splitting into odd and even 
parts, we arrive at the following equalities, to be used subsequently: 

vtoA = 7rQX±H and Avro = —H^A^vro. (5) 


5 



2.2 An alternative notation for the unattenuated case 

We now introduce some notation that emphasizes the L^(M)-duality arising in the Pestov- 
Uhlmann reconstruction formulas [29]. The main novelty below is the introduction of I±. The 
general unattenuated transform can be defined over functions / G L?‘{SM) as 

f{^pt{x,v))dt, (x,v)Gd+SM. (6) 

Jo 

Considering integrands of the form / G Lp‘{M) (i.e. /(x, u) = /(x)), and X_\_h for some h G 
Hq{M), let us define the unattenuated transforms 

Iof-=If, I±h:=I{X^h), 


(that is, in the definition of Iq, we identify / with its pullback by the canonical projection 
SM —)• M). These transforms are continuous in the following settings when {M,g) is non¬ 
trapping 

Jo : L\M) ^ Llid+SM), : H^{M) ^ Ll{d+SM), 

where is a weighted space with weight ^(x, u) = g{ux,v). Define the scattering relation 
a : dSM —)• dSM as follows: if (x, u) G d+SM, a{x,v) = G d-SM, and if 

(x, u) G d-SM, a{x,v) = ¥?_t-(x,-u)(x, u) G d+SM. Recall the definitions of A± : L‘j^{d+SM) —)• 
L^^|(d5M) and their adjoints A^, introduced in [29]: 


A±w{x, 9) 


w{x,v) on d+SM, 
±woa{x,v) on d-SM, 


A*^u := {u±uo a)\d+sM- 


The fundamental theorem of calculus along a geodesic reads 


IXu = {no a - u)\q_^_sm =-A*-U. (7) 

With V’ : SM —)■ d-SM the endpoint map 'ip{x,v) = (x, u) and h G L^{d+SM), we denote 

hy := h o a o Ip the function extended by free geodesic transport to SM, i.e. solution of the 
equation 

Xu = 0 {SM), u\3_^sm = h. 

Straightforward computations using Santalo’s formula yield that for h G L‘j^{d+SM), 

I^h = 2TT{h^)o, I^h = -2Tr{X_ih^)o. ( 8 ) 
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The V± decomposition. Let us define the antipodal scattering relation ai : d^SM —?■ d^SM 
as {v I—)• —v) o a, and write L‘j^{d+SM) as the direct sum V+ © V_, where h G V+ (resp. V_) 
iff h is even (resp. odd) with respect to the involution ai. Since a function of x only can be 
regarded as an even function of v on 5M, and a vector field can be regarded as an odd function 
of V on 5M, it is straighforward to establish that 


Range Iq © and Range I_\_ C V_. 

Moreover, we have the following lemma. 

Lemma 2.1. The direct sum L?^{dj^SM) = V+ © V_ is orthogonal. 

Proof. For h G L^(d+5M), define u = where r = /o(l) (1 denotes the constant function 

equal to 1 on M). Santalo’s formula allows to show that the map h i—>■ u is L‘j^{d+SM) —)• Lp‘{SM) 
continuous and ||u|| 2 , 2 ( 5 ^) = \\h\\i,2^(^Q^sM)- Moreover, u is even/odd in v whenever h G V+/V_, 

so that, if h G V+ and g G V_, and upon calling u = ^ and v = have 

{h, g) = {u,v)l 2 i^sm) = 0, 

hence the proof. □ 


Inversion of Iq and I±. We now revisit the inversion of the operators Iq and I±, previously 
established in [29], adapted here to the present notation. Recall the notation (/ G L‘^{SM)) 
for the solution of a transport problem of the form 

Xu = -f (SM), u\ 9 _sm = 0, (9) 

and for / G define Wf = {X±U'f)o. It is established in [29] that W extends as a 

smoothing (hence compact) operator W : -G and that the L^-adjoint of W is 

given by W*h = 

Proposition 2.2. Let {M,g) a simple surface with boundary. Then we have for every f G 
L‘^{M) and every h G Hq{M), 

f + W^f = ^Ilw, w = ^A\H.A.Iof, (10) 

h + {W*fh = -^IqW, w = ^AlH+A-I^h. (11) 

ZTT 4 

Remark 2.3. Formulas (10) and (11) differ slightly from [29, Theorem 5.4] because it is stated 
there that the solution to the transport problem 

Xu = -/, u\a+sM = w, 

is + w,p, which is what the Fredholm equation there is based upon. The correct answer would 
be u^ + {w — lof)^, which in turn yields the modified formula (10). 
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Proof. Proof of (10). Let / G C°°{M) and define as in (9) so that u\g^sM = lof- Applying 
itqX = itoX±H (derived in (5)) to (9), we obtain 

/ = vro/ = -TToXu^ = -TToX±Huf = -ttoX±HuL = -{X±HuL)o, 

so / can be obtained from the last equation if we can relate Hu_ to the known data lof- In 
order to do so, we use the commutator [H,X] to write a transport equation for Hui_ 

XuL = -f => X{HuL) = (Axu^)o = -IP/, 


so that Hu[_ satishes the transport problem 


X{HuL) = -Wf, HuUa.SM = ^H.A.Iof\9_sM ■■= ri, 


which means that 


Hu[_ = + w^, w = T] o a. 


Upon applying (Ai_L-)o, we obtain 

{X^Hui)o = W^f-^Ilw. 

Since we have established that / = —{X±HuL)o, we conclude that 


/ + w = Q(iL_^_/o/)|a_SM^ o a- 

A* —A* 

In terms of operators, we can also write w as w = — -H-A-Iof. Moreover, it can 

be seen that the function A^H-A-Iq/ has V+ symmetry, so it is annihilated by Thus, 
Equation (10) follows, extended to every / G L?‘{M) by density of C°°{M) in Lp‘{M). 

Proof of (11). Let h G C°°{M) with h\QM = 0, and let solve the transport equation 


Xu = —Xj_h {SM), 'u\g_sM = 0. 

Upon projecting onto odd functions of v, we have Xu^^^ = —X±h. Direct manipulations and 
the use of the commutator formula imply 

Xu^^^ = -X^h X{Hu^^’^ -h) = -X^W*h -= -X^W*h. 

Moreover, the trace of is given by u^-^^\gsM = ^A-I±h. Since the function 

satishes the transport problem 

X{Hu^^^ -h) = -X^W*h, Hu^^%_sm = ^H+A_I^h\g_sM, 



we deduce that 


— h = ^ + w^, where w := (^Hj^A-I_\_h\Q_sM^ ° cn. 

Upon averaging over 0 and rearranging terms, we obtain 

h + (W*fh = -^I^w. 

A* —A* 

In terms of ^4^ operators, w can also be written as w = 5 ^ 2 — ~H^A-Ij_h. We can see that 
A*_H^A-I±f has V_ symmetry and as such is annihilated by /q. Thus, Equation (11) follows, 
extended to every h G Hq{M) by density. □ 

Remark 2.4. Inspection of symmetries shows that A-I^f in (10) is odd in v and A-I±f in 
(11) is even in v. This tells us that the expressions H^A^Iof and H^A-I±f are redundant, as 
one eould just replace H± by H. This further emphasizes the similarity between formulas (10) 
and (11), for both of which the operator Aj^HA- acts as a first step in the postproeessing of 
data, though on a different subspaee depending on the formula. 

2.3 Construction of explicit continuous right-inverses for Jq and If 

The question of surjectivity of If and If have proved to be crucial for answering boundary rigid¬ 
ity questions (see [30] where a proof of surjectivity of If appears) and constructing holomorphic 
integrating factors in a prior study of the attenuated transform (see [34] where Lemma 4.5 states 
that If is surjective), which are so important to the present approach. Such proofs relied on 
pseudodifferential arguments on an extended simple compact manifold, and did not construct 
explicit preimages of either operator. In order to derive and implement explicit inversions, con¬ 
structing explicit preimages becomes a necessity, and we notice here that, while writing formulas 
(10) and (11) in a way that emphasizes duality, we also notice that the right-hand sides involve 
If and If directly. Under the assumption that the operators Id + W‘^ and Id + (VU*)^ are 
invertible, this allows for an explicit construction of continuous right-inverses of If and If. 

Remark 2.5. Although it would be enough to show that Id + lU^ is injective, which is open at 
present for general simple surfaces, it is shown in [13] that the operator W admits a bound of 
the form ||lU||j;^ 2 ^j ^2 < C||Vk||oo- This implies that if curvature is close enough to eonstant, the 
operators Id + IT^ and Id -T (W*)^ are invertible via Neumann series. Whether this qualitative 
assumption covers the case of all simple surfaces remains open at present. 

Proposition 2.6. Suppose the operators Id + and Id + are lf{M) —)• L‘^{M) in¬ 

vertible. Then for every A: G N, the operators iij_ : C^{M) —^ C^{d+SM) and Rq : C^+^(M) —)• 
C^{d+SM) defined by 

R± := ^AfH.A.Ioild + W^)-\ Ro := -^A\H+A_I^{Id + (12) 
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are continuous and satisfy I^R±f = f and IqRoJ = f for f smooth enough. 

Proof. Suppose that Id + W'^ and its adjoint are invertible. Then their inverses map any C^{M) 
to itself. This is because the kernels of Vk, IT* are proved in [29] to be smooth so that, e.g., if 
/ G L?‘{M) solves / + IT^/ = /i, where /i G C^{M) and W‘^ f is smooth, then f = fi — W‘^ f 
is . Since the operators W, W* are also L? —)• continuous, then similar arguments allow to 
show that {Id + and its adjoint are C^{M) —)• C^{M) continuous. 

The relations I'fRi, = Id and IqRq = Id are straightforward to check, as a directly applica¬ 
tion of equations (10) and (11) and the invertibility of Id + W'^ and Id + (kk*)^. 

It remains to prove that 

\\R±f\\c'={d+SM) < C\\f\\ck(^M)^ \\^og\\ci^{d+SM) < C\\g\\ck+ii^M)- (13) 

Looking at the compound expression of these operators, we see that A- and ^4^ preserve 
norms since the scattering map is smooth and the function r(x, v) is smooth on d^SM whenever 
dM is strictly convex (see [35, Lemma 4.1.1 p.ll5]), H± preserve norms as convolution 
operators, and Iq : C^{M) —)• C^{d+SM) and /j_ : C^+^(M) —)■ C^{d+SM) are continuous since 
the geodesic flow is smooth. □ 

Remark 2.7. A study of symmetries with respect to the involution ai shows that p = R±f and 
q = Rog thus constructed satisfy p G V_ and q ^ V+. This is compliant with the continuity 
statements (13), as any component of p in V+ would be annihilated by and any component 
of q in V- would he annihilated by Jq . 

2.4 Holomorphic solutions to certain transport eqnations 

A crucial tool in the inversion of attenuated ray transforms is the construction of holomorphic 
integrating factors, whose existence relies on the surjectivity of Iq and I{[. In the simple Rie- 
mannian setting, it is proved in [26, Theorem 4.1] that the transport equation Xu = —F (for 
F G C°°(S'M)) admits holomorphic solutions if and only if F is of the form F = fi + X±f 2 
for some functions /i ,/2 G C°°{M). Although uniqueness of such solutions may not hold (e.g. 
adding a constant to such a solution makes another solution), a constructive approach, inspired 
in part by [26, Theorem 4.1], is to look for an ansatz, holomorphic by construction, of the form 

u = {Id + iH)p.^jj + {Id + ill)q.^, 
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where p is a smooth element in V_ and g is a smooth element in V+, so that is odd and 
is even. Plugging this ansatz into Xu = —fi — X±f 2 , we obtain 

Xu = X{Id + iH)p^ + X{Id + iH)q^ 

= -i[H, X]p^ - i[H, X]q^ {Xp^ = Xq^ = 0) 

= -i{X±p^)o - iXjAp^- iX±{q^)o 

Therefore, a sufficient condition for u to solve Xu = —F is if p and q satisfy 

= and =/a, 

which we may solve explicitly using the previous section. Using Proposition 2.6, we summarize 
this construction in the following result, whose proof is straightforward and omitted. 

Proposition 2.8. Under the hypotheses of Proposition 2.6, and given A: G N, /i G C^{M) and 
/2 G the function 

u = 27ri[{Id + iH){R^h)^ - (Id + iH){Rof 2 U, 

is a holomorphic solution of Xu = —/i — X±f 2 on SM, satisfying the estimate 

\W\\c’^{SM) < C'(ll/l|lc'=(M) + ll/2|lc''+i(M))) 

where Rq and R± are defined in (12). 

3 Inversion of the attenuated ray transform over functions and 
vector fields d la Kazantzev-Bukhgeim 

In [12], the authors provide reconstruction formulas for functions and vector fields from knowl¬ 
edge of their ray transforms in the case where the metric is Euclidean and the domain is the 
unit disk. The present section generalizes these ideas to the case of simple Riemannian surfaces. 

3.1 Reconstruction of a function 

In this first approach, we follow the idea in [12] that, if we can find a solution u*, holomorphic 
with Mq = 0 of Xu* + au* = —/, then projecting this equation onto Hq gives 

/ = -V+U-i - V-ut - auQ = 
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after which one must explain how to express Ui in terms of known data. Here and below, 
assuming that the operators Id + and Id + {W*)‘^ are invertible, we denote Iq^± the ray 
transform restricted to integrands of the form /i + X±f 2 , where /i ,/2 £ C°^{M). If D : = 
lo/i + I±f 2 = I{fi + X±f 2 ), we can reconstruct /i ,/2 (and in turn, /i + X±f 2 ) from D by 
carrying out the following steps: 

1. Compute the V+/V- decomposition of D. This decomposition will be given by h 
a^D) £ V±, with ai the antipodal scattering relation. 

2. Reconstruct /i from the projection of D onto V+ by inverting (10). 

3. Reconstruct /2 from the projection of D onto V- by inverting (11). 

In what follows, we summarize the above procedure by writing Iq]^D = /i + ^_l/ 2 - 

Theorem 3.1. Let {M,g) a simple Riemannian surface with boundary and f,a £ 

Then f is uniquely determined by its attenuated geodesic transform via the reconstruction formula 

f = 2ir?_(Im(e-h:^))i, h' := - w')\a^sM, (14) 

where w is an odd, holomorphic solution of Xw = —a, D = {{Id — iH){e~'^{u\gsM))- with 
''J‘\d-SM = 0 and u\a^sM = laf, o.nd w' is a holomorphic solution of Xw' = —lQ]_{AfD), given 
by Proposition 2.8. 

Proof. Step 1: find n* a holomorphic solution of Xu* + au* = —f. Call u the solution 
to Xu + au = —/, u\q_sm = 0 so that u\q^sm = laf ■ Let w a holomorphic, odd solution of 
Xw = —a. Then v = e~''’u solves the transport problem 

Xv = e~'''{Xu-{Xw)u) = e~'''{Xu +au) =-e~''’f, w|a+SM = e""'Ia/, v\ 9 _sm = 0, 

(15) 

The right-hand side b := e~'^f is holomorphic. This is because, since the product of holomorphic 
functions is holomorphic, so are (convergent) powers series of holomorphic functions. Moreover, 
plugging the expansion w = wi + w^ + ... {wj £ Hj) into the exponential, we see that the 
fiberwise decomposition of 6(x, v) looks like 

b = f{l + w + ^ + ...) = + 

bo bi 

02 

SO bo = f. We now want to decompose v into v = v* + h'^, where v* is holomorphic and h'^ is 
constant along geodesics, then we will have that v* solves Xv* = —b. We proceed as follows. 
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First decompose v = ^), where = {Id ± iH)v. The function ^ solves the 

transport equation 

Xv^-^ = X{v - iHv) 

= {Id-iH)Xv + i[H,X]v 
= —{Id — iH)b + i{X_\_v)o + iX^VQ 
= -(/ - + iX±_VQ. 

Integrating along geodesics, we deduce that 

h{f — i{^±'v)o) — iI±vo = A^v^~'^\dsM, 

where the right-hand-side A*_{v^~'^\dsM) = A*_D, where D has the expression in the statement 
of the theorem and is known from data laf ■ Therefore the right-hand side {f — i{X±^v)Q)—iX^VQ 
can be reconstructed upon inverting Jq and /_l, a relation which we denote 

(/ - i{X^v)o) - iX^VQ = I^X^*_D. 

Let w' a second holomorphic function such that Xw' = —(/ — i(^±u)o) -k iX±vo, constructed 
following Proposition 2.8, that is, w' = 2'jri{Id + iH){p^ -|- q^) with p, q solving 

I1p= {f-i{X±v)Q) and = iuo- 

With w' thus constructed, we have X{v^~'> —w') = 0, i.e. {v^~^ —w') is constant along geodesics. 
Upon rewriting v as v = + w') + ^{v^~'^ — w'), we see that the first term is holomorphic 

and the second is constant along geodesics. In other words, v is of the form v = v* + h'^, where 
V* = 5 ( 1 ’^"''^ + w') and h'^ = ^{v^~'^ — w'). Additionally, with this choice of w', we have 

^0 = + w'o) = + 27rf(g^)o) = ^(^^0 + ilU) = 0- 

Finally, defining u* = u*e"', we see that u* is holomorphic as the product of holomorphic 
functions, and, using the last equation, we see that Ug = (u*e’^)o = Ug = 0 , and that it solves 

Xu* + au* = -be^ = -/. 

Projecting the equation above into flo) we obtain 

-/ = d+uli + r?-< + «< ^ 

so it remains to show how to compute u\ in terms of known data. 

Step 2: express in terms of known data. We now write 

= u{ — = Ui — {u*)i = 2 z(Im(u*))i = 2 z(Im(n* — u))i, 
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where the first equality comes from the fact that u* ^ = 0 and the last comes from the fact that 
u is real valued. We now write, by definition, 

u* — u = (v* — v)e^ = 

so that Im(M* — u) = —Im(e"'/i^), and we arrive at the expression 

u* = —2i(Im(e"'/i^))i. 

Now looking to compute h', since we proved that h'^ = — w'), then 

h' = - w')\d^SM 

= ]^v^~\d+SM - ^(27ri)[(Ici + fF)(p^ + g^)]|a_^5M, 

where r^p = g^ I*q = ivQ, hg - iI±vo = A*_{v^~^\dSM)- 

The proof is complete. □ 

Remark 3.2. This proof generalizes the one completed in [12] in the Euclidean case when the 
domain is a disk. In order to eomplete the argument there, it is required to relate the fiberwise 
Hilbert transform with the Hilbert transform of the domain for the so-called divergent beam 
transform (see [12, Lemma 4-1]), whieh in turn uses the singular value decomposition (SVD) of 
that operator, established in earlier references (see [12[ for detail and references there). 

While this SVD, specifie to the choice of metric and domain, does not seem straightforward 
to generalize systematically, the present eonstruction of holomorphic solutions allows here to 
simplify the proof by bypassing these steps altogether. 

In isothermal coordinates (x, 0), Equation (14) takes the following form; 

/ = -g_ut = -e-^^d (< e") , <(x) = ^ 6)) dO, (17) 

with d := \{dx + idy). Constructing both holomorphic solutions w,w' using Proposition 2.8, we 
arrive at the following reconstruction procedure: 

1. Construct w = 27ri(/ + iR)(R_i_a)^, odd holomorphic solution of Xw = —a. 

2. Compute v\^_^_sM = e^Hf and v\a_sM = 0, then v‘'~1\qsm = {Id - iH){v\dsM)- 

3. Reconstruct g and vq from log — iI±vo = A'*_{v^~'l\gsM)- (inversion of Iq and I±) 

4. Construct p and q from il[_p = g and I^q = ivo- (inversion of Iq and If, see Sec. 2.4) 

5. Construct h' from (16). 

6. Reconstruct u* then / according to (17). 
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3.2 Reconstruction of vector fields 

As in [12], the method of proof above can easily be generalized to the case of reconstruction of 
vector fields, which we now present. Such a problem has applications to Doppler tomography 
in media with variable index of refraction, previously studied in [11, 12], Integrands of the form 
/(x, 9) = /o(x)+/i(x) cos(0—a) are also considered in [2] in the Euclidean setting. In our context, 
a smooth real-valued vector field takes the form V(x, u) = {Fi +F-i)/2, with F±i an element of 
and Fi = E_i. In isothermal coordinates, we may write Ei(x, 9) = (/i(x) — i/ 2 (x))e“'^^’^)e*® 
for some real-valued functions /i,/ 2 , and this would correspond to integrating the vector field 
/i(x)da; f2{^)dy. 

Theorem 3.3. Let {M,g) a simple Riemannian surface with boundary and V(x, u) = (Fi -|- 
F_i)/2 a smooth vector field as above with F_i = Fi. Then at every x G M where a(x) fi 0, 
F_i (and hence V) can be uniquely reconstructed from data laV via the formula 

F-i = -4ir/_ Qr/_(Im(e"'/i[^))i^ , h' := - w')\d^sM, (18) 

where w is an odd, holomorphic solution of Xw = —a, D = {{Id — iH){e~'^{u\QSM))-)\dSM with 
'F\d+SM = laV and u\q_sm = 0? and w' is a holomorphic solution of Xw' = 

Proof. Start from the equation 

Xu + au =—{Fi + F_i)/2, u\q^sm = laV, u\q_sm = 0. 

Step 1: find u* a holomorphic solution of Xu* + au* = —V. Let w = wi + + ... be 

an odd, holomorphic solution of Xw = —a and define v = e~'"u. Then v satisfies the transport 
problem 

Xv = -e-"'(Fi + F_i)/2 = -G(x, v) = -(G_i + Go + Gi + G 2 ...), 

where G_i(x,u) = F_i(x,u)/2 and Go(x) = —wiF-i/2. Split v into an holomorphic and 
antiholomorphic part, i.e. look at u = where = {IdF H)v. v^~'i satisfies the 

transport equation 

Xv^-^ = X{v - iHv) = {Id - iH)Xv + i[H, X]v 

= -{Id - iH)G + i{X^v)o + iX^VQ 
= —(Go — i(Aj_u)o) + iX±vo — G_i. 

The data D defined in the statement is D = v^~'^\dsM- Using the Hodge decomposition, we now 
write G_i = Xg + X±h for g, h two functions on M, where g fulfills the additional prescription 
g\dM = 0. Then the previous transport equation becomes 

X{v^~'> +g) = -(Go - i(A_Lu)o) -h X±{ivo - h). 
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Note that {v^ '> + g)\d_^_sM = ^\d+SM = (Id - iH){e '^Iaf)\d+SM is known, so that upon 
integrating the transport equation along each geodesic we can see that the data gives us 

/o(Go - i(X^u)o) - I±(iuo -h) = AtD, 

with D defined as in the statement of the theorem. Now construct w' a holomorphic solution of 
Xw' = -(Go - i{X^v)o)+X^{ivo -h) = 

so that + g — w') = 0, i.e. there exists h defined on d^SM such that + g — w' = h^. 

We now decompose v = + w' — g) + — w' + g) = v* + h'^, where the first 

term is holomorphic and the second is constant along geodesics. In particular, we get that 
Xv* = —G{x,v) where v* is now holomorphic. Then defining u* = we find that u* is 

holomorphic and satisfies 

Xu* + au* = -(T_i + Fi)/2. 

Looking at the projections onto and flo yields the relations 

g-UQ = -F_il2, g-Ui + auQ = 0 , 

which implies the reconstruction formula, at each point where a does not vanish: 

F_i = 2 t/_ . 

Step 2: obtain from the measurements. This part is, again, similar to the proof of 
Theorem 3.1. We write 

u\ = u{ — = Ui — (u*)i = 2z(Im(u*))i = 2z(Im(n* — u))i, 

where the hrst equality comes from the fact that u* ^ = 0 and the last comes from the fact that 
u is real valued. Next we have the relation 

u* — u = (v* — v)e^ = —e^hX 

so it remains to compute which after unrolling definitions, 

d' = h'Xa+SM = -w' + g)\d+SM = “ w')\d_^_sM, 

where both terms are, again, compatible from data: v^~'>\q^sm = D\d+SM and w' is a holomor¬ 
phic solution of Xw' = —Iq^[A*_D], a solution of which can be explicitely constructed following 
Proposition 2.8. This ends the proof. □ 
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4 A second approximate formula for functions, conditionally in¬ 
vertible via Neumann series 

While Theorem 3.1 reconstructs functions exactly and, in some sense, in a “one-shot” fashion, 
it presents a couple of weaknesses: (i) it involves all values of rc(x, 9) throughout SM, which 
in turn involves storing three dimensions of data when everything should be dealt with using 
two-dimensional structures, and (ii) the effects of curvature (which need iterative corrections as 
in [18]) are not being corrected in the right place. 

We now propose an algorithm based on a more direct interplay of the Hilbert transform with 
transport equations as in [19, 29], which leads to a Neumann-series based inversion, faster in 
implementation, and valid when curvature is close enough to constant and the attenuation is 
small enough in norm. Such an algorithm is then implemented in Section 5. 

We first state a result about a certain family of operators generalizing the operator Wf = 
{X±U'^)o first defined in [29]. These operators appear as error operators of the next reconstruc¬ 
tion formula. We relegate the proof and some remarks about these operators to the appendix. 

Proposition 4.1. Let h G L°^{SM) such that Vh G C^{SM). Then the operator Kh : L?‘{M) —)> 
L‘^{M) defined by K^f := (X_i_n-^^)o is well-defined and continuous, and there exists a constant 
C such that 


WKhh^^L^ < CdlVKllooljhlloo + ll^hllci). (19) 

We now state the main result of the section. 

Theorem 4.2 (Iterative inversion). // / G Lp‘{M) and laf denotes its attenuated GXRT with 
given attenuation a G C‘^{M), then the function f satisfies the following equation 

f + Kf = ^Ilp, rj=^AlH{e-^U)_, ( 20 ) 

where {e~^Iaf)- denotes extension of e~^Iaf from d+SM to dSM by oddness w.r.t. v i-;- —v, 
w is an odd, holomorphic function solving Xw = —a and K : L?‘{M) —>• L?‘{M) is a bounded 
operator satisfying the following estimate 

< C (IIVk]]^ + ell“ll-(llalloo||VKlloo + jjalb)) . (21) 

Proof. Proof of (20). Let w a holomorphic, odd solution to Xw = —a. If u is the solution 
to Xu + au = —f with boundary condition u\q_sm = 0, then the function v = e~'^u solves 
the transport problem Xv = —b := —fe~'^ with boundary condition v\q_sm = 0) so that v is 
no other than ™. Applying vroA = tt(}X±H (derived in (5)) to the equation Xv = —b, we 
obtain 


-/ = -vro6 = ttqAu = {X^Hv)o = {X^Hv.)o. 
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The task now is to make Hv- more explicit. Hitting the transport equation satisfied by v with 
the Hilbert transform H, we write a transport equation for Hv 

X{Hv) = -Hb - X^vo - {X±v)o. 

As b is holomorphic, we have Hb = fH{e~^) = — 1) and upon taking the even part of 

the last equation w.r.t. 0 i—)■ 0 + vr, the function Hv- solves the transport equation 


X{Hv-) = if{cos\iw — 1) — (Ai_i_n)o, 
where we have used that w{yi, 0 + vr) = —rc(x, 6). This tells us that 

Hv- = ^ ^{X^v)o ^ 


where is constant along geodesics. We can deduce r] from the fact that the first two terms in 
the last r.h.s. vanish on d-SM, so that r/ = {Hv-\d_sM) ° ot, and since V- is known from data 
on dSM, so is rj. More precisely, we have, for any (x, v) G dSM 


v-{x,v) = -(e '^Iaf)-{^,v) := 


lg-«,(x,.)j^J(x,„) 

^e-“’(X’-")/a/(x, -v) 


if (x, v) G d^SM, 
if (x, v) G d-SM. 


( 22 ) 


Upon applying the operator vroA^, we obtain 


f = -{X^Hv-)o 

= i))o - WiX^v)o - {X^Vi.)o 

= i))o - 

i.e. this equation takes the form 

f + Kf =-{X±r]^)o, rj = {Hv-)\o_sm ° a, 
where Kf := W{X^u^^~")o - 


Note that upon rewriting ™ ™ the operator K can be rewritten as 

Kf = W{X^u^)o + 

= W^f + 

Equation (23) now corresponds to (20) upon noticing that Ifr] = —27r{X±r].^)Q, and that 


lU = II ( \{H{e-'^IJ)-)\a_sM oa]=Il Qfl+^(F(e-4/)_) 
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where the last equality comes from the fact that A*_ applied to an odd function on dSM makes 
a function with V+ symmetry, which in turn is annihilated by 

Proof of (21). The theorem will be proved once we show that the operator K satisfies the 
bound (21). From the last equation, we see that K = + WK^-^ — iK^^, where hi = e~'^ — 1 

and /i 2 = cosh re — 1 and the notation Kh refers to Proposition 4.1. It is proved in [13] that 
||IF||l 2^^2 < CIIVkIIoo for some constant C. By virtue of Proposition 4.1, we deduce that for 
i = 1,2 

< C'(||VK||oo||hj||oo + lll^^illci)- 

Now we bound ||hj||cxD < and ||P/ij|| 0 i < for 

j = 1,2, and together with the fact that w is constructed following Prop. 2.8, it satisfies 
estimates of the form 


||u;||cfc < C'ljallcfe, A: = 0,1, 2,... 

Combining all these estimates together, we obtain estimate (21). □ 

Remark 4.3. In the absence of attenuation a = 0, equation (20) is exactly the Fredholm equation 
(10), in which case wfK,v) = 0 so that Kf = W'^f and the right-hand side of (23) is the post¬ 
processing of unattenuated data. 

Corollary 4.4. Under the hypotheses of Theorem ^.2, if ||Vk||oo and ||a||(;2 are small enough 
so that estimate (21) implies ||iir|| 2 , 2_,.^2 < 1, then f can he reconstructed from equation (20) via 
the Neumann series 

oo 

n=0 

Remark 4.5. As in the unattenuated case, this restriction on ||Vk||oo and ||a||(;2 is of rather 
qualitative nature and does not tell us whether all simple cases will work, and whether this 
approach would work for attenuations less than C^. This is to he contrasted with successful 
numerical reconstructions below, which work for both discontinuous attenuations and cases of 
metrics arbitrarily close to non-simple. 

5 Numerical implementation 

We now present a brief implementation of an inversion of (20) via a Neumann series approach. 
We use the code developed by the author in [18] for the unattenuated case, augmented with 
attenuation. The domain M is the unit disk {x = (x, y),x‘^ + < 1} endowed with the metric 

g(x) = where 

A(x) = 0.2 ^exp ^ ^ ^ 
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describing a region of “low sound speed” near xq and “high sound speed” near —xq. The effect 
on geodesic curves can be seen Fig. 1. For such a domain, it can be computed that the boundary 
is strictly convex and there are no conjugate points. 

The influx boundary d+SM = x (—f, f) is parameterized by x(/l) = and ingoing 

speed direction 6{(3,a) = P + n + a (in this case, /? + vr is the direction of the unit inner normal). 
M is represented into the unit square [—1,1]^, discretized in an equispaced cartesian fashion 
with N = 300 points along each dimension. 



Figure 1: Geodesics cast from the boundary point x(/3) = with, from left to right, 

/? = 0, ^, |, ^. The colorbar on the right indicates the values of the background sound speed 
c(x) = 


The function / and the attenuation a are displayed on Fig. 2. Note that both quantities 
contain jump singularities. 



Figure 2: From left to right: Unknown /, attenuation a, function n defined on d^SM such that 
I^n = a (so that w = 2'Ki{Id + iH)n^ is a holomorphic solution of Xw = —a). Axes for n are 
{I3,a) G [0,27r] x (^, f). 


Strategy for inversion. Equation (20) is of the form 


f + Kf = Lalaf, 


(24) 
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where laf represents forward data, LaD = is an approximate inversion 

and we assume that K is a contraction. Once la and La are discretized (call their discretized 
versions with the same name for simplicity), discretizing K separately and computing a finite 
sum of to reconstruct / may introduce additional numerical errors due to the 

fact that (24) may not be satisfied at the numerical level. A better approach is to set directly 
—K = Id — Lala and to implement a finite sum of the series 

OO 

/ = J^(/d - LalafLalaf. 

fc =0 

We now briefly explain how the operators la and La are implemented. 


Forward operator la- The computation of the forward data is done by discretizing the in¬ 
flux boundary d+SM into an equispaced family {IIi,aj)i<i< 2 N,i<j<N and for each data point, 
we compute Iaf{/3i,aj) by discretizing the system of ODEs over t G [0,Tjj] where Tij = 


X = I 


i/cos0\ 0 = g '''W(d A(x) COS0 — 93 ;A(x) sin0), 
Vsiny/ 

with initial conditions x(0) = ^(0) = Pi + tt + aj, along which we compute 

lafiPuOij) = / /(x^,,a^.(t))exp ( / a{x 0 ^^aj{s)) ds] dt, 


via a discrete sum, with the solution of the ODE above. 


Computation of a holomorphic solution w of Xw = —a. Following Proposition 2.8, we 
look for w in the form w = 2Tri{Id + iH)n^, where n solves I^n = a. This requires implementing 
n = R±a, with R± = ^A^H-A-Io{Id + . Following [18], {Id + is computed via 

a few iteration of a rapidly convergent Neumann series, Iq is the unattenuated X-ray transform. 
In the present case, n is represented on the right-hand plot of Figure 2. Once n is computed, as 
the expression of La only involves values of w on dSM, then we can compute 

w\dSM = {27ri{Id + iII)n^)\0SM = 2.-Ki{Id + iII){n^)\osM = 2m{Id zi?)A+re, 

where the Hilbert transform is processed via Fast Fourier Transform on the columns of A+n. 
As n has V_ symmetry, A+re amounts to extending n to d-SM by oddness w.r.t. 6 ^ 0 + tt 
(the latter is much more straightforward numerically). 
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Approximate inverse L^. For D a data function defined on a discretization of d^SM, we 
wish to compute LaD = The computation of w and H is explained in the 

previous paragraph. A^h{f3, a) is computed by combining values of h{(3, a) and h at the endpoint 
of the geodesic starting from coordinate (/3,a). The main technical step is the computation of 
/^, which in isothermal coordinates can be simplified as follows (see [18, Sec. 3.1.2]): 

where Vx- is just divergence on a cartesian grid. Integrals in can then be discretized using 
finite sums and the divergence is implemented using finite differences. 

Both computations of and require computing several geodesic endpoints, which is the 
main bottleneck of the code. An alternative option, trading memory for much shorter CPU time, 
is to compute and store all endpoints required at first, and reusing them in further Neumann 
iterations. 



Figure 3: Left: attenuated data laf (data for Experiment 1) with /, a from Fig. 2. Right: 
unattenuated data Iq/ for comparison. The differences may be observed in magnitude but not 
in support. 


Experiments. We present two experiments, in which a and / refer to the functions displayed 
on Eigure 2. 5a refers to the function M 9 x i—)• 5a (x). 

• Experiment 1. (low attenuation) Neumann series based reconstruction of / from /„/. 

• Experiment 2. (high attenuation) Neumann series based reconstruction of / from I^af- 

Experiment 1 successfully and stably reconstructs / within 3 Neumann iterations, as shown in 
Eig. 4 (up to numerical inaccuracies, and given the fact that jumps can never be fully captured 
exactly). As this example works even if a is discontinuous, this is to be contrasted with the 
regularity requirements on a from Theorem 4.2. 
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Figure 4: From left to right: reconstruction frc from /„/ and pointwise error \ frc — f\ after one 
iteration, frc from /«/ and pointwise error \ frc — f\ after three iterations. 


Experiment 2 displays a divergent Neumann series, due to the fact that attenuation 5a is too 
high. This is in agreement with the smallness requirements of Corollary 4.4 on the attenuation 
coefficient, as the operator norm of the error operator in (21) potentially grows like 
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Figure 5: From left to right: I^af (data for Experiment 2), pointwise error \frc — f\ after one 
iteration, pointwise error \frc — f\ after two iterations. Divergence of the algorithm ensues due 
to the appearing unstable modes. 
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A A certain family of operators - proof of Proposition 4.1 

On simple surfaces, it is proved in [29] that the operator Wf = {X±uf)Q is smoothing on simple 
surfaces, so that the equation reconstructing a function from its unattenuated ray transform is 
of Fredholm type. 

It has been observed that if / is now a function on SM, the corresponding operator no 
longer has such properties. However, the operators Kf^ introduced in Proposition 4.1 appear 
naturally as error operators in Theorem 4.2, and they generalize W since W = Kh with h = 1. 
In general, may no longer be smoothing, but we can still obtain Lp‘{M) —)• Lp‘{M) continuity 
with estimates on the norm in terms of h and the ambient curvature. We must recall some facts 
about Jacobi fields (or variations of the exponential map), following [16]. For ^ G we 

may decompose along the frame {X{t), X±{t),V{t)} at the basepoint tpt{yi,v) as 

d(ptiC) = Ci(x,u,t)X(t) + C2(x,u,t)X_L(t) +C3(x,u,t)P(t). 


The structure equations then provide us a differential system in t for the coefficients Cg 

Cl =0, C 2 + Cs = 0, Cs — K(7x,i;(t))C2 = 0. 


In particular, we may express the variation fields dipt{X±) = aX±{t) — dV{t) and d^ptiV) = 
—bX±(t) + bV{t) in term of two functions a, b defined on P, solving 


d + K(7x,t;(t))a = b + = 0, 


a b 
a h 


(o) = [S?] 


The assumption of simplicity implies that b does not vanish outside {t = 0} and is thus positive 
for all t > 0. Note the constancy of the Wronskian ab — bd = 1. 


Proof of Proposition 4-1 ■ For a smooth function 4>{x,v) on SM vanishing on dSM, we first 
write, using the chain rule, 

^ t ( x , h ) pT{x,v) 

x± / dt = / {a{x,v,t)X±(j){(ft{x,v)) - d{x,v,t)V4>{(pt{x,v))) dt. 

Jo Jo 

We then rewrite (keeping (x, u) implicit), for t ^ 0: 


ab 


1 , 


aX±4> — dV(j) = -{bX±(j) — bV4>) + —d + — V(f = — - V{(f o ipf) -\- - V(f. 


Integrating this equality for t G (0,r(x, u)) and integrating by parts over 5x (using that 4> 
vanishes at dSM), we arrive at the conclusion that 




(j){(pt{x,v)) dt dS{v) 
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Now replacing (j){x,v) by f{x)h{x,v) and using the fact that V{fh) = fV{h), we arrive at 

Kh{f){y^) = ^ 

= Kh,if{^) +Kh,2f{^), 


upon expanding the sum. The operator K^ i is just as well-behaved as the operator W and for 
the same reason: defining q{x,v,t) := t) ^ ( Ux’v’t) ) ’ shown in [13] that |q'(x, u,t)| < 
CIIVkIIoo for every {x,v,t) £ V. So we can rewrite 

Kh,if{x) = -^f [ q{x,v,t)h{ipt{x,v))f{'y^^y{t)) b{x,v,t) dt dS{v) 

J Syi J 0 

= / 9(x,w(y),t(y))/i((^t(y)(x,u(y))/(y) dMy, 

so that the kernel kh,i{^, y) = (?(x, u(y), t(y))/i((^j(y)(x, u(y)) of Kh,i is bounded (hence in L‘^{M x 
M)), i.e. the operator K^ i : Lp‘{M) —)■ Lp‘{M) is continuous (in fact, compact) with an operator 
norm controlled by C'||VK||oo||/i||oo- On to the study of the second term 

^h,2/(x) = / ,2/ ’ ^^ /(7x,^.(t)) b dt dS{v), q2{x,v,t) := Vh{(pt{x, v)). 

^'^Js^Jo b^{x,v,t) 

The function q 2 satisfies Jg q2(x,v,0) dS{v) = fg Vh(x,v) dS{v) = 0, so that the integral is 
expected to make sense as a principal value integral. Note that near t = 0, we have b{x,v,t) = 
t + t^c{x, V, t) where c is smooth on V, and 1 + t^c{x, v, t) does not vanish on D since b does not 
vanish outside {t = 0} by simplicity of the surface. More precisely, we write 

‘^T^Kh, 2 f{^) = [ [ b dt dS{v) = Af{x) + Bf{x) + Cf{x) 

Js^Jo onx,u,t) 

where, upon writing b{x, v,t) = t + t^c{x, v, t), we define 


^/(x) 

5/(x) 

C/(x) 


'Sx Jo 


{^ 1 ^} Qryiy^ y Qj 

y ’ fh.At))bdtdSiv) 


Is^ Jo 


r(x,t>) g2(x,U,t) - q2{x,V,0) 

i2- - 


/(7x,«(i)) b dt dS{v) 


f ^ ^,c{x,v,t ){2 + t^c{x,v,t)) u\\ i. ^or \ 

/ / 92(x,u,t - f-, , .2 ( -Ivl2- b dt dS{v). 

is^Jo {i + t^c[x,v,t)y 


Upon changing variable {v,t) i—>■ y{v,t) = 7x,u(i) (with Jacobian dMy = b dt dS{v)), the C term 
becomes an operator with bounded kernel, i.e. —)• bounded, with operator norm controlled 


25 









by Ik 2 ||oo) be. ||Tb/i||oo- On to the B term, we may write \q 2 {^,v,t) — g 2 (x,r;,0)| < 
uniformly on 2?, so that the kernel of B has an integrable singularity and the B term becomes 
an operator with integrable kernel, i.e. —)• bounded, with operator norm controlled by 

||y/i||ci. On to the A term, we assume without loss of generality to be working in isothermal 
coordinates. We change variable {v,t) i—y(u,t) = 7x,'u(i) to make appear 


Am = [ 

JM 


g2(x,u(y),0) 


/(y) iM, = f dy, 

Jm idg[^,y)r 


Im K(x,y))2 JM y^g\ 

where dy now represents the Lebesgue measure on M^. Expansions near x give that 

dg(x,y) = e^(y^|x-y| + |x - ypdi(x, y), 

^ lx_y| + 

this allows to rewrite ^4 as a Calderon-Zygmund operator of the form 


bl/(x) 



g2(x,g,o) 

|x-yp 


/(y) dy + 



g3(x, y) 
|x-y| 


/(y) dy, 


9 


x-y 

x-y|’ 


(25) 


where is uniformly bounded by CIlE/iUci. By virtue of [17, Theorem XI.3.1], the first term 
together with the zero mean value condition Jgig 2 (x, 0,0) d0 = 0 is an operator —)• 

continuous, with an operator norm bounded by Csup^g^ || 92 (x, •)IIl2(§i)) in turn bounded by 
C||iy^||oo- The second term of (25) is another weakly singular operator whose operator norm 
can be bounded by C'HE/iUjji as well. □ 
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